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Collective states of interacting non-Abelian anyons have recently been studied mostly in the context of certain 
fractional quantum Hall states, such as the Moore-Read state proposed to describe the physics of the quantum 
Hall plateau at filling fraction = 5/2. In this manuscript, we further expand this line of research and present 
non-unitary generalizations of interacting anyon models. In particular, we introduce the notion of Yang-Lee 
anyons, discuss their relation to the so-called 'Gaffnian' quantum Hall wave function, and describe an elemen- 
tary model for their interactions. A one-dimensional version of this model - a non-unitary generalization of the 
original golden chain model - can be fully understood in terms of an exact algebraic solution and numerical 
diagonalization. We discuss the gapless theories of these chain models for general su(2)fc anyonic theories and 
their Galois conjugates. We further introduce and solve a one-dimensional version of the Levin- Wen model for 
non-unitary Yang-Lee anyons. 
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I. INTRODUCTION 

Over thirty years ago Leinaas and Myrheim pointed 
that in systems confined to two spatial dimensions parti- 
cles with exotic exchange statistics, more general than those 
of bosons and fermions, are possible. Such particles with 
arbitrary exchange statistics were later coined anyons by 
WilzeclP. Today it is widely believed that this possibility is 
indeed realized in the fractional quantum Hall effect. An even 
more intiiguing form of statistics has recently received con- 
siderable attention, namely that of non-Abelian statistics, first 
proposed in a seminal paper by Moore and Read^. This form 
of statistics can occur in two-dimensional systems in which 
introducing excitations gives rise to a macroscopic degener- 
acy of states. Upon braiding the excitations, the wave func- 
tion (or better, the vector of wave functions) describing the 
system, does not merely acquire an overall phase, but can ac- 
tually transform into one another, as described by a unitary 
braid matrix acting within the degenerate manifold. In gen- 
eral, these braid matrices do not commute, hence the name 
non-Abelian statistics. While several systems have been theo- 
retically proposed to exhibit quasiparticles with non-Abelian 
statistics, such as unconventional + ipy superconductors"*, 
rotating Bose-Einstein condensates'', or certain heterostruc- 
tures involving a novel class of materials, so-called topologi- 
cal insulators*, the system currently being cast under intense 
experimental scrutiny is the fractional quantum Hall effect ob- 
served at filling fraction v ~ 5/2, with some evidence sug- 
gesting that this state is indeed non-Abelian in nature- 

An early attempt to describe this v = h/2 quantum Hall 
state came in the form of the so-called Haldane-Rezayi wave 
function', which can be thought of as a d-wave paired BCS 
condensate of composite fermions - electrons with, in this 
case, two flux quanta attached. A peculiar feature of the 
Haldane-Rezayi state is that its gapless edge modes are de- 
scribed by a non-unitary conformal field theorjESG^I with cen- 
tral charge c = —2. However, non-unitary dynamics can- 



not describe a physical system as it would violate basic piin- 
ciples of quantum mechanics. It has therefore been argued 
that in general the non-unitary nature of an edge state indi- 
cates that the underlying phase is not bulk gapped, but in fact 
critical'^ in which case the edge states lose their identity as 
they dissipate into the gapless bulk. For the Haldane-Rezayi 
state it indeed turns out that it does not describe a gapped topo- 
logical phase, but rather the gapless state at a quantum phase 
transition between two gapped states - the Moore-Read quan- 
tum Hall statd^ and a so-called strong pairing quantum Hall 
phase^. Since the work of Haldane and Rezayi many other 
wave functions have been proposed, which appear to have 
non-unitary edge state theorie !r^"^ i One particularly well- 
known example is the so-called 'Gaffnian state These pro- 
posed states typically arise from numeiical studies of model 
Hamiltonians, for which it is oftentimes hard to determine 
whether they describe a gapped or gapless phase. The ultimate 
fate of these newly proposed states has therefore remained an 
active field of research. 

In this paper, we will study the excitations of the putative 
quantum liquid described by such non-unitary wave functions 
from an 'anyon-model' perspective. Thus, we explore the 
physics that would occur if one were to assume that such 
excitations indeed exist and study in detail their collective 
behavior that would result. We refer to these excitations as 
non-unitary anyons and model the interaction between these 
anyons in the same way as it was done in the unitary casJ^H. 
It will turn out that the Hamiltonians one is led to consider for 
the non-unitary anyons are, not surprisingly, non-Hermitian. 
Nevertheless, the energy spectrum turns out to be completely 
real. One can hope that the results of such a study would shed 
some light on the questions surrounding the physical mean- 
ing of non-unitary anyons. Of course, in the context of two- 
dimensional classical statistical mechanics models having a 
non-unitary conformal field theory (CFT) desciibing a critical 
phase is a situation which often occurs - see, for instance, the 
(restricted) height models described in Ref. [2j which are gen- 
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eralizations of the 'restricted solid-on-solid' modelJ^U- and is 
not problematic. 

Before closing this introduction with the outline of the pa- 
per, we should mention one additional point, which moti- 
vates the study on non-unitary anyon models: the concept of 
'Galois-conjugation' of conformal field theories. In short, the 
action of the Galois group relates different CFT's which have 
the same number of fields, obeying the same fusion rules. The 
non-unitary models we encounter in this paper, are in fact 
Galois-conjugates of unitary CFT's, and thus, the study on 
non-unitary anyons is closely related to the study of Galois 
conjugation in conformal field theorjBSEll Por the purposes 
of this paper we will however not need to use any details of 
these connections. 

This paper is organized as follows: In Sec.|ll]we discuss the 
Yang-Lee anyonic chains. We start by giving a review of the 
derivation of the Hamiltonian of the (ordin ary, i .e. unitary) 



golden chain^^ of Fibonacci anyons in Sec. 



11 A 



and a short 



mathematical description of the Galois-conjugated model, the 



so-called Yang-Lee chain in Sec. II B including an explicit 



derivation of their non-Hermitian, microscopic Hamiltonians, 
and discuss their relation to the 'golden chain' models of 
Ref. [23]via Galois conjugation. We continue in Sec. II C 
with a discussion of the algebraic structure underlying these 
microscopic Hamiltonians, which in turn allows for an analyt- 
ical identification of their gapless theories in terms of certain 
non-unitary minimal models of conformal field theory. We 



next present a number of exact, numerical results in Sec. II D 



and discuss the topological symmetry protecting the critical 



states in Sec. II E before rounding off with a discussion of the 



case of general level k in su(2)fe theories in Sec. II F 



In Sec. Ill we turn to 'doubled Yang-Lee' models and, in 
particular, a non-unitary generalization of the (unitary) high- 
genus ladder model which was previously studied for Fi- 
bonacci anyons in Ref. '28^, where we introduce the model in 
Sec. |III A| We then discuss the phase diagram in Sec. [Ill B| 
present the analytical solution at special critical points in Sec. 
I III CI and finish with exact numerical results in Sec. IIII Dl 

The final section summarizes the results and discusses their 
relevance to hypothetical non-unitary topological phases. An 
appendix contains a detailed discussion of the conformal en- 
ergy spectra at the four critical points discussed in the two 
models. 



II. YANG-LEE CHAINS 

At the focus of this manuscript are anyonic models that 
are certain non-unitary generalizations of unitary non-Abelian 
anyon models, which have been extensively studied in the re- 
cent pasPsElHlIl The basic constituents of the generalizations 
considered here are non-unitary, non-Abelian anyons. Like 
their unitary counterparts they carry a quantum number that 
corresponds to a generalized angular momentum in so-called 
su(2)fc anyonic theories, which are certain deformationJ23 of 
SU(2). We first concentrate on an elementary example where 
there is only a single anyon type by explicitly considering the 
anyon theory su(2)3. In the unitary version of this theory the 



elementary degrees of freedom are often referred to as 'Fi- 
bonacci anyons', and it is their non-unitary counterparts which 
we term 'Yang-Lee anyons'. We will return to a discussion 
of the non-unitary generalizations for general su(2)fc anyonic 
theories in Section HTFl 

We start by quickly reviewing the basic construction of mi- 
croscopic (chain) models of interacting non-Abelian anyons, 
following the ideas of the 'golden chain' model of Ref.|23]and 
the detailed exposition of Ref. 3 1 The construction of these 
models proceeds in two steps. First, we describe the general 
structure of the Hilbert space of these models in a particular 
'fusion chain' representation, which is identical for the uni- 
tary and non-unitary models. In a second step we turn to the 
microscopic form of the Hamiltonian capturing interactions 
between the anyons. While this second step is quite similar 
for the unitary and non-unitary cases, the microscopic Hamil- 
tonians for the two cases are distinct. 



A. The golden chain 

The elementary degrees of freedom in our microscopic 
model are the particle types (or generalized angular momenta) 
of the su(2)3 anyonic theory. In its simplest form (considering 
only integer momenta) this theory contains a trivial particle 
(or vacuum state), which we denote by 1, and an anyonic par- 
ticle, which we label as r. These particles can form combined 
states according to the fusion rules 



1x1 = 1 



1 X T = T 



T X T = 1 



(1) 



The non-Abelian nature of the anyonic r-particle reveals it- 
self in the multiple fusion outcomes when combining two of 
these particles. Our chain model then consists of L such t- 
particles in a one-dimensional arrangement as depicted on the 
top of Fig. [T] where L denotes the number of sites of the 
chain. Since pairs of r-particles can be fused into more than 
one state, such a system of L non-Abelian anyons spans a 
macroscopic manifold of states, i.e. a vector space whose di- 
mension grows exponentially in the number of anyons. It is 
this manifold of states that constitutes the Hilbert space of our 
microscopic model. To enumerate the states in the latter we 
define a so-called 'fusion chain' as illustrated on the bottom 
of Fig.[T] Here the original r-particles constituting the chain 
are denoted by the lines which are 'incoming' from above. 
The links in the fusion chain carry labels {xi} which again 
correspond to the particle types of the su(2)3 theory. Reading 
the labels from left to right a labeling is called admissible if 

T T T T T 



xq Xi xi xs 



FIG. 1. A chain of Fibonacci or Yang-Lee anyons (denoted by the 
t's in the upper row). The set of admissible labelings {xi} along 
the fusion chain (lines) constitutes the Hilbert space of the Yang-Lee 
(and Fibonacci) chains. 



at each vertex the fusion rules ([TJ of su(2)3 are obeyed, i.e. 
a r label is followed by either a 1 or t label, while a 1 la- 
bel is always followed by a -r label. Every such admissible 
labeling then constitutes one state in the Hilbert space of our 
anyonic chain. Considering periodic boundary conditions, i.e. 
XL — it is straight forward to show that the dimension of 
the Hilbert space is given in terms of Fibonacci numbers as 

dimL = FibL_i + Fibi+i , 

where Fibi denotes the i-th Fibonacci number, defined by 
Fibi+i = Fibi + Fibi_i and the initial conditions Fibi — 
Fiba = 1. 

We now proceed to the second step of our construction, 
the derivation of a microscopic Hamiltonian. In doing so we 
follow the perspective of the original 'golden chain' model^^ 
in assuming that interactions between a pair of neighboring 
T particles - mediated, for instance, by topological charge 
tunnelinJ^H- will result in an energy splitting of the two possi- 
ble fusion outcomes in Eq. ([T]). Our Hamiltonian captures this 
splitting by projecting the fusion outcome of two neighboring 
T particles onto the trivial fusion channel, i.e. assigning an 
energy of — —1 to the fusion of two t particles into the 
trivial channel and an energy of Et = to the fusion into the 
r channel. This anyonic Hamiltonian is thus reminiscent of 
the common Heisenberg Hamiltonian for SU(2) spins, which, 
for instance, projects two ordinary spin-l/2's onto the singlet 
channel and assigns a higher energy to the alternative triplet 
channel. 

To explicitly derive the Hamiltonian in the Hilbert space of 
fusion chain labelings introduced above, we note that in this 
basis the fusion of two neighboring t particles is not explicit. 
To get direct access to this fusion channel of two neighbor- 
ing T particles, we need to locally transform the basis as de- 
picted in Fig. |2] The matrix describing this transformation 
is typically called the F-symbol, which can be thought of as 
an anyonic generalization of Wigner's 6j-symbol for ordinary 
SU{2) spins. Its general form (in the absence of fusion mul- 
tiplicities) is given in Fig. [3] 

Assuming that we know the explicit form of the i^-symbols 
(see the next section for more details), we can now explicitly 
derive the microscopic Hamiltonian in the fusion chain basis. 
After the basis transformation, the fusion channel of the two 
neighboring anyons is manifest, so by means of a simple pro- 
jection we can assign an energy to each of the fusion channels. 
The final step left after this projection, is to transform back to 
the original basis, which again employs the F-symbol. 

To make the individual steps of this derivation more ex- 
plicit, we consider the example of Fig.|2]in more detail. Let 
us specify the five possible labelings of three neighboring fu- 
sion chain labels Xi-i, Xi, Xi+i, where in Fig. |2]we depicted 

T T T T T T T T T T 

I I I I I . F I I I 

X\ Xi X-i . . . Xq X\ X-i . . . 

FIG. 2. The F-symbol describing the local change of basis. 
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FIG. 3. The general form of the F-symbol. 
the case where the site label is i = 2, 

1) , |1,x,t) , \t,t,1) , |t, , \t,t,t)} . 

After performing the basis transformation shown in Fig.|2] the 
following labels satisfy the fusion rules at each vertex and thus 
form the new basis 

{|1, 1, 1) , |1, T, t) , |t, T, 1) , |t, 1, t) , |X, T, t)} , 

where Xi is the fusion channel of the two neighboring t par- 
ticles. In the transformed basis, we can project onto the trivial 
channel, by means of a projection Pi x, where the subscript i 
denotes that we are acting on anyons i and i + 1, while the 
label 1 denotes we are projecting onto the 1 channel. So, the 
part of the Hamiltonian acting on anyons i and i + 1, which 
we denote by acts on the Hilbert space as 

W \xi^i,Xi,Xi+i) = 

(2) 

Here, we have used that for the su(2)fc anyonic theories we are 
considering, the F-symbols are their own inverses. Moreover, 
we projected onto the 1 channel, which we favored, because 
of the overall minus sign. The total Hamiltonian then simply 
becomes the sum of (|2]i over all positions 

L 

H = Y,H\ (3) 

i=l 

where we assume periodic boundary conditions, i.e. xl = xq. 

To describe the Hamiltonian of the various types of anyon 
chains we consider in this paper, we only have to specify the 
explicit form of the F-symbols (apart from the fusion rules, 
which determine the Hilbert space). The explicit form of the 
Hamiltonian then follows from Equation (|2|. 

B. Galois conjugation and non-unitary models 

Now that we have expressed the Hamiltonian in terms of the 
F-symbols, we should explain how to obtain the F-symbols 
for a given anyon theory. As stated, the F-symbols transform 
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FIG. 4. The g-deformation parameters of Fibonacci and Yang-Lee 
anyons correspond to different primitive roots of unity. 



between two different fusion bases as illustrated in Fig. [3] 
As such, the exact form of the these symbols can be deter- 
mined self-consistently by identifying a circular sequence of 
basis transformations, which yield a set of strongly overcon- 
strained nonlinear equations called the 'pentagon equations' 
(for a more detailed exposition see, for instance, Refs.[3T|and 
[38T l. While finding a solution to these pentagon equations is in 
general a highly non-trivial task, it has been shown that they 
allow only for a finite set of inequivalent solutions, a prop- 
erty which goes under the name of 'Ocneanu rigidity', see for 
instance Ref. 39 For the su(2)fe anyonic theories of interest 
here, the complete set of possible F-symbols can be found, 
e.g., in Ref. 00] where they were obtained by using quantum 
group techniques. 

The different i^-symbols are found to have a general form 
that depends on a single, so-called 'deformation parame- 
ter' q only. This deformation parameter has to be chosen 
appropriate 1)0^ and it turns out that for the su(2)fc anyonic the- 
ory it must be one of the [k + 2)"^^ primitive roots of unity, 
i.e. of the form 



,p-27ri/(fc+2) 



(4) 



where the integer index p runs from l<p<(A: + 2)/2 (and 
p and fc + 2 are relative prime). The process of increasing 
the index p by one, i.e. going from one root of unity to the 
next, is what is usually referred to as Galois conjugation. For 
our example theory, su(2)3, we can thus identify two possible 
values for q, which are illustrated in Fig. |4] These two Ga- 
lois conjugated theories corresponding to deformation param- 
eters, q = e^'^'/^ and q = e^^^/^, then precisely correspond 
to the cases of Fibonacci and Yang-Lee anyons, respectively. 
The explicit form of the i^-symbols and in particular the non- 
diagonal 2x2 matrix for can then be writterP^in terms 
of this deformation parameter q as 



(5) 



For Fibonacci anyons we set q = e^^*/"'', in which case g ^ + 
1 + q = I + 2cos(27r/5) = (1 + Vb)l2 = is the golden 
ratio, and the i^-symbol becomes the unitary matrix 



^Fib 



-1/2 



6-1/2 



(6) 



The golden ratio, of course, is one solution of the equation 
= 1 + x, which is an algebraic analog of the fusion rule t x 
X = 1 + T of the su(2)3 anyonic theory. The process of taking 
the Galois conjugate of the original Fibonacci anyon model 
corresponds then simply to the substitution — > — 1 /(/), where 
— !/(/) is the other solution to the equation — 1 +x. In terms 
of the deformation parameter q, this amounts to choosing the 
other possible value of q — e^'^*/^, which indeed yields q^'^ + 
1 + gi = 1 + 2 cos(47r/5) = -1/0. The F-symbol for Yang- 
Lee anyons thus becomes the (invertible) non-unitary matrix 



Yang— Lcc 



-Z(/.l/2 



1/2 



(7) 



Having obtained the F-symbols in both the unitary as 
well as the non-unitary case, we can now write down the 
Hamiltonians for the Fibonacci and Yang-Lee chains. On 
the states \xi^i,Xi,Xi+i) e {|1, r, 1) , |1, -r, r) , |x, r, 1)} 
both Hamiltonians act in the same diagonal way, W = 
diag{— 1, 0, 0}. Acting on the states \xi^i, Xi, Xi+i) G 
1, r) , |t, t, t)|, the Hamiltonians take the following 
formgSl 



TTl 

Fibonacci 



H 



Yans— Lc 



0-3/2 ^-1 
02 J03/2 
103/2 _^ 



(8) 



Before discussing these anyonic models in further detail, 
we note that while Galois conjugation changes some aspects 
of these models, i.e. the parameters in their respective Hamil- 
tonians get 'Galois conjugated', this turns out to be a rather 
mild change, since the underlying algebraic structure of these 
models remains largely untouched. As a consequence, the 
non-unitary Yang-Lee chains allow for an analytic solution 
similar to their unitary counterparts as first obtained for the 
'golden chain' model in Ref. .23, We will discuss the details 
of this analytical solution and its resulting gapless theories in 
the next Section. 



C. Algebraic structure and analytical solution 

In the original golden chain paper^^*, it was shown that the 
Hamiltonian ([3]) based on the unitary F-symbols (|6]), can be 
exactly solved. Following a similar sequence of steps as in the 
unitary case we show here that all non-unitary models allow 
for an exact analytical solution as well. We briefly outline 
these steps in the following, for a more detailed discussion 
see the references 12313 II 

As a first step, it was noted that the operators i.e. the 
summands of the Hamiltonian, form (upon suitable normal- 
ization) a known representation**^ of the Temperley-Lieb al- 
gebra'*-' with "d-isotopy"- parameter d, namely = —dH^, 
where d — (f). These operators satisfy the Temperley-Lieb 
algeba 



d ei 
for 



(9) 
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FIG. 5. The allowed 'height' configurations of the chain models. 
This diagram is the Dynkin diagram of the Lie algebra Ak — su{k + 
1). The Fibonacci and Yang-Lee anyons correspond to fc = 3. 



To show that this is the case, we will first map the Fibonacci 
and Yang-Lee chains onto so-called restricted solid-on-solid 
(RSOS) or 'height' models^^. To do this, we first note that 
one can relate the Fibonacci anyons to the anyonic theory 
su(2)3, a theory with four particles (see section 11 F for the 



more general case of su(2)fe), which can be labeled by their 
'spin' j = 0, 1/2, 1, 3/2. The fusion rules of these particles 
are given in table |l] One immediately notices that the parti- 
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TABLE I. Fusion rules of the su(2)3 theory. 

cle with I ~ 1 has the same fusion rules as the Fibonacci 
anyon. Indeed, the Fibonacci anyon model can be viewed as 
the 'integer subset' of the su(2)3 theory. In addition, one no- 
tices that fusing a particle j with the particle 3/2, one finds 
j X 3/2 = 3/2 — j, swapping integer 'spin' to half integer 
'spin'. This allows one to map the Fibonacci chain, consisting 
of 'spin'-l particles, to a chain of 'spin' -1/2 particles, by fus- 
ing the labels X2i-i of the odd sites with the particle 3/2. In 
this way, the new labels are constant labelings of a chain con- 
sisting of 'spin'-l/2 particles. Performing this map is advanta- 
geous, because the transformed Fibonacci-chain can now di- 
rectly be mapped onto a height model, in which the allowed 
heights take the values 0, 1/2, 1, 3/2, which can be seen as the 
nodes of the Dynkin-diagram A4, which is given in Fig. [5] 

The well-known Jones representation of the Temperley- 
Lieb algebra, acting on the Hilbert space of the (transformed) 
Fibonacci chain, can now be obtained by using the so-called 
modular S-matrix."*^ Explicitly, this representation is known 
to take on the following form: 



Y *-'0,a;i_i Y ^O.XiJ^i 



(10) 



where Sj^^j^ is the modular 5* matrix of the su(2)3 theory, 
and the labels ji , j2 are the 'spins' of the corresponding parti- 
cles, which in general take the values 0, 1/2, 1, . . . , fc/2, and 



in terms of these, the S'-matrix elements read 

■(2ji 



fc + 2 



1)(2J2 



1)tt 



(11) 



In equation ( [TOj i, the row-label ji of the 5-matrix elements 
which appears here takes on the value ji = 0, which is the cor- 
rect value in case of the (unitary) Fibonacci chain. In taking 
the Galois conjugate, one has to replace this row-label ji = 
by the label ji = 1 (corresponding to a r particle [see also 
section fll F| ). 

The action of the e^'s on the height states, precisely cor- 
responds to the action of the part of the Hamiltonian acting 
on the ket \xi-i,Xi, Xi+i), which shows that the operators of 
Hamiltonian of the Fibonacci chain do indeed form a repre- 
sentation of the Temperley-Lieb algebra. 

We will now employ this observation, and map the Fi- 
bonacci chain to an integrable, two-dimensional statistical 
mechanics model, namely an RSOS model, which in the 
present case is based on the Dynkin diagram A4 - see Fig. 
[5] In an RSOS model, the degrees of freedom are the heights 
(or nodes of the Dynkin diagram). These heights live on the 
vertices of the lattice, with the constraint that heights of neigh- 
boring vertices correspond to nodes of the graph which are 
linked. 

The two-row transfer matrix of the RSOS-model (depicted 
in figure|6| can be written in terms of the plaquette weights of 
the square lattice, namely T = T2T1, with 



Ti=]JW[2i] T2=]JW[2z + l] 



where the plaquette weights are of the form 



(12) 



sin(^) ^-^sin(J^) 



(13) 



(14) 



With if = (n,<5.^,., 



the identity operator. 




FIG. 6. The two-row transfer matrix of the RSOS models. 



The new feature of the non-unitary models consists in the 
fact that equations ( [T3| l and ( fl?] ) for the plaquette weight in 
the transfer matrix now involve a general integer p = 2ji + 1 
labeling the primitive root of unity d?]) beyond the unitary case 
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of p = 1 (corresponding to ji = 0) studied earlier For a 
general value of p, the Hamiltonian associated with this lattice 
model can be obtained by taking the extreme anisotropic limit, 
namely u — > 0+ (see for instance 1441). One has. 



T = cxp(-a(H + ci) + 0(a2)) 



(15) 



with a = 



^ 1 and ci an unimportant constant. 



which indeed gives that the Hamiltonian is of the form H = 
— gX^i'^i' establishing the mapping from the chain to the 
RSOS models. 

We can now use the known results about the phases of 
the RSOS models, to find the behavior of the chain models. 
In the case of the Fibonacci anyons (i.e., when fc = 3 and 
p — 1, which appear in the plaquette weights W and the pa- 
rameter a), the anti-ferromagnetic golden chain Hamiltonian 
can be obtained by taking the limit u with u positive. 
The corresponding (exactly integrable) RSOS model^^ is crit- 
ical and known^^ to be in the universality class of the tri- 
critical Ising model described by the minimal model 7W(4, 5) 
of conformal field theory. For the ferromagnetic golden chain, 
which is obtained from the RSOS models with < \ u\ ^ 1 
and u negative, the critical behavior then turns ouP^lII to be 
in the universality class of the parafermions corresponding 
to the minimal model variant"*^ A^(5, 6). 

In the case of the Galois-conjugated Yang-Lee anyonic 
chains (corresponding to fc = 3 and p = 2), the same se- 
quence of steps as above results in a mapping to another fam- 
ily of exactly integrable RSOS models^ and in this cas e the 
gapless theories turn out be non-unitary minimal modelJ^^El 
In this non-unitary case, the Hamiltonian can be obtained from 
the two-row transfer matrices, but in this case, both the anti- 
ferromagnetic as well as the ferromagnetic chains are obtained 
in the limit u — > with positive u. The difference in the 
Hamiltonians stems from the sign-changes in the isotopy pa- 
rameter d, namely from positive (but not necessarily bigger 
than one) in the anti-ferromagnetic case, to negative in the 
ferromagnetic case. The critical behavior in these two cases is 
described by the (non-unitary) minimal models A^(3, 5) and 
M{2,5) respectively. 



D. Numerical results 

We have numerically studied the excitation spectra of the 
Yang-Lee chains by exact diagonalization of systems with up 
to L = 32 anyons, typically using periodic boundary condi- 
tions. These excitation spectra not only allow for an indepen- 
dent identification of the conformal field theory describing the 
gapless collective state, as discussed in the previous section, 
but also reveal further details about the correspondence be- 
tween continuous fields and microscopic observables. In par- 
ticular, the low-energy states of a conformally invariant sys- 
tem can be identified with conformal fields and the excitation 
spectrum is expected to take the form 
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FIG. 7. Conformal excitation spectrum of the 'antiferromagnetic' 
Yang-Lee chain. The spectrum matches the non-unitary minimal 
model M(3, 5) with central charge c = —3/5, which is often re- 
ferred to as 'Gaffnian' theory. Primary fields /, o, e, ip of this con- 
formal field theory are indicated by squares, descendant fields by 
circles. We also indicate the 'topological flux' of each energy eigen- 
state, which indicates the topological symmetry sector. 



where h and h are the (holomorphic and anti-holomorphic) 
conformal weights of a given CFT with central charge c. Ei 
is a non-universal number, v a non-universal scale factor, and 
L the length of the chain. To match the excitation spectra of 
the Yang-Lee chains to these CFT predictions we consider the 
family of so-called minimal models M{p, p') (where p and p' 
are mutually prime) with central charge 



1 



6{p_~p') 
pp' 



/^2 



and conformal weights 

h{r, s) = 



[rp — sp')^ ^ {p ^ p'Y 
App' 



(17) 



where the indices r and s are limited to 1 < r < p' and 
1 < s < p. We note that the labels (r, s) and {p' — r,p ~ s) 
coiTespond to the same field. 

In the following, we will discuss our numerically obtained 
excitation spectra for 'antiferromagnetic' and 'ferromagnetic' 
couplings, which are plotted in Figs. |7] and [8] respectively. 

The antiferromagnetic chain.- We first turn to the 'anti- 
ferromagnetic' chain, for which the pairwise anyon-anyon in- 
teraction energetically favors the trivial fusion channel 

r X r — > 1 . 

The conformal field theory describing the critical behavior of 
this model is the non-unitary minimal model A^(3,5) with 
central charge c = —3/5, which is also referred to as the 
'Gaffnian' theory'^. The four primary fields of this CFT and 
their respective scaling dimensions A ^ h + haie 
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with the non-trivial fusion rules 



12 



X e : 
e X e 



(7 

I 



a X ip = e 
e X — a 
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For completeness, we give the conformal dimensions of the 
fields (with minimal model labeling) in table and note that 
this model is a particular Galois conjugate of the su(2)3 CFT. 

To identify the gapless theory numerically, we typically 
perform the following procedure: We first look at the two 
lowest energy eigenvalues in the spectrum, Eq and Ei, and 
by identifying the energy gap AE = Ei — Eq with the dif- 
ference of the two lowest scaling dimensions we can identify 
the non-universal scale factor 2ttv /L in ( [T6| ), which we subse- 
quently set to 1 thereby rescaling the entire energy spectrum. 
This identification of the two lowest energy eigenvalues with 
conformal operators also allows to identify an overall energy 
shift, e.g. setting the energy of the trivial operator / with scal- 
ing dimension h{l, 1) + h{l, 1) = to zero. In the case at 
hand, there is only one negative scaling dimension, so the low- 
est energy corresponds to — 2/imin = —1/10, while the second 
lowest state corresponds to the identity operator, with zero en- 
ergy. At this point, all the energies are fixed, and indeed the 
rescaled and shifted numerical spectrum is found to reproduce 
the position of the (other) primary fields (indicated by green 
squares in FigjTll, as well as the descendants (indicated by red 
circles in Fig. [tJi. 

The ferromagnetic chain- We now turn to the 'ferromag- 
netic' chain, for which the pairwise anyon-anyon interaction 
energetically favors the r-fusion channel 

T X T ^ T . 

The critical theory is the non-unitary minimal model M{2, 5) 
with central charge c = —2 2/5, which is commonly referred 
to as the 'Yang-Lee' theorj^^H. The two primary fields of 
this CFT and their respective scaling dimensions A = h + h 
are 



(19) 
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with the non-trivial fusion rule 



e X e 



I + e 



Again we give the conformal dimensions of the fields (with 
minimal model labeling) in table|ll]for completeness. 



M{2,5) 



h{r,s) 


s = 1 


2 


h{r,s) 


s = 1 


r = l 





3/4 


r = l 





3 


1/5 


-1/20 


3 


-1/5 



TABLE n. Kac table of conformal weights for the non-unitary mini- 
mal model M{3,5) and (2, 5) . We only displayed fields with odd 
r labels, to avoid duplicates. 
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FIG. 8. Conformal excitation spectrum of the 'ferromagnetic' Yang- 
Lee chain. The spectrum matches the non-unitary minimal model 
M(2,5) with central charge c = —22/5, which is commonly re- 
ferred to as 'Yang-Lee' theory. The primary fields I, e of this con- 
formal field theory are indicated by squares, descendant fields by 
circles. We also indicate the 'topological flux' of each energy eigen- 
state, which indicates the topological symmetry sector. 



The spectrum of this theory, after the appropriate shift and 
rescaling of the energy is displayed in figure [8] which beauti- 
fully reproduces the primary fields, as well as descendants to 
a high level, constituting the spectrum of the Yang-Lee model. 



E. Topological symmetry 

Before considering further generalizations of the Yang- 
Lee chains we mention another peculiarity of these anyonic 
chains. Like their unitary counterparts the Yang-Lee chains 
exhibit an unusual, non-local symmetry. This symmetry, 
which was dubbed a 'topological symmetry' in the context 
of the golden chain modeP^l corresponds to the operation of 
commuting a t particle through all particles of the chain. The 
so-defined operator, which we denote by Y , is found to com- 
mute with the Hamiltonian and for the su(2)3 theory has two 
distinct eigenvalues, thus defining two symmetry sectors. Its 
matrix form is given by 



L-l 



''0' ■ 



(20) 

This definition solely in terms of the F-symbols immediately 
suggests a generalization of this symmetry to the case of the 
non-unitary Yang-Lee models studied above by simply replac- 
ing the F-symbols with their non-unitary counterparts (|7|. For 
both the unitary and non-unitary variants the two eigenvalues 
of the respective topological symmetry operator are yi — (p 
and 2/2 = —^l<t>- In the unitary case these are identified as 
no-flux / r-flux symmetry sectors. This assignment is simply 
reversed in the Galois conjugated, non-unitary case. 

For the unitary models it has been demonstrated that this 
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topological symmetry protects the gapless ground state of the 
interacting anyon chain model^'*; it was shown that all rel- 
evant operators (in a renormalization group sense) that have 
otherwise identical quantum numbers as the ground state, e.g. 
the same momentum, fall into different topological symme- 
try sectors. We have performed a similar symmetry analysis 
for the Yang-Lee chains at hand. For chains with either anti- 
ferromagnetic or ferromagnetic couplings, we have evaluated 
the eigenvalue of the topological symmetry Y for all eigen- 
vectors of the Hamiltonian and thereby assigned topological 
symmetry sectors to the primary fields of the conformal field 
theory describing their energy spectrum. These assignments 
are given in Figs. [7] and [8] for antiferromagnetic and ferromag- 
netic chain couplings, respectively. A situation similar to the 
unitary models emerges: For both signs of the coupling the 
ground state with conformal dimension h + h<Ois found to 
be in the topologically non-trivial (or r-flux) sector, while all 
other primary fields with the same momentum are found to be 
in the topologically trivial (or no-flux) sector. 

For the unitary anyon chains this topological protection 
mechanism has subsequently been cast in a broader physi- 
cal picture^"* interpreting the gapless modes of the anyonic 
chains as edge states at the spatial interface of two topolog- 
ical liquids, and the conclusion that anyon-anyon interactions 
result in a splitting of the topological degeneracy for a set of 
non-Abelian anyons and the nucleation of distinct topologi- 
cal liquid within the parent liquid of which the anyons are 
quasiparticle excitationJSEl The similarity of our results for 
the topological symmetry properties of the non-unitary anyon 
chains thus raises the question whether a similar interpretation 
would also hold for the non-unitary systems at hand. 

F. su(2)fc theories and generalized Galois conjugation 

While we have focused our discussion in the preceding 
sections on the Fibonacci / Yang-Lee anyonic theories for 
which fc = 3, our results can be readily generalized to a 
much larger class of anyonic theories, so called su(2)fc Chern- 
Simons theories. These theories, which are certain quantum 
deformations^* of SU(2), have a finite number of representa- 
tions which are identified with the (anyonic) particle content 
of the theory. For a given (integer) deformation parameter 
k these representations/particles are labeled by a generalized 
angular momentum 

. _ 1 k-1 k 

The fusion rules for combining two of these momenta (or their 
associated particles) then resemble the tensor product rules of 
ordinary SU(2) angular momenta, and explicitly read 

min( ji +j2 -j-2) 
jl X j2 = XI • (21) 

33=\'ji-h\ 

For all fc < 2 the first non-trivial fusion rule (i.e. a fusion rule 
with multiple fusion outcomes) is encountered for the particle 



carrying generalized angular momentum j = 1/2, for which 
one finds 

^xi=0 + l. (22) 

Upon further inspection of the fusion rules (see Table |l] for 
an example) one notes that for anyonic theories with odd k 
one can map all half-integer representations onto integer rep- 
resentations by fusing the half-integer representations with the 
highest representation of the theory, the so-called 'simple cur- 
rent' with angular momentum |. In particular, this mapping 
identifies the fusion rules of the particle carrying angular mo- 
mentum j = 1/2 with those of the particle carrying angular 
momentum j = {k — l)/2 and Eq. ( |22] i becomes 



Note that for odd k the fusion rules are now given entirely 
in terms of integer representation^^. For the su(2)fc=3 theory 
this mapping thus turns the fusion rules of the j = 1/2 particle 
into those of the j — 1 particle, i.e. 1x1 = + 1, which 
is the momentum representation of the non-trivial fusion rule 
in Eq. ([T]). One can thus equally identify the x particle of 
the Fibonacci / Yang-Lee theories with both the j ~ 1 and 
j — 1/2 angular momenta in su(2)3. 

The su(2)fe generalizations of the Fibonacci / Yang-Lee 
chains, which we want to consider in the following, are chains 
of particles carrying angular momentum j = 1 /2, which we 
can alternatively identify with particles carrying angular mo- 
mentum j = {k — 1) /2. The microscopic chain model is then 
constructed in an identical way to what we outlined for the 
case of the Fibonacci / Yang-Lee theories in Section [ll] i.e. 
we first build a Hilbert space by considering admissible label- 
ings of the fusion chain basis for an incoming j = 1/2 (or 
j = {k — l)/2) particle and then derive a Hamiltonian via 
the local projection onto one of the two fusion outcomes in 
Eqs. ( |22] i or ( |23] l, respectively. Our notion for the sign of the 
interaction also remains unchanged, with 'antiferromagnetic' 
interactions favoring fusion of neighboring particles into a to- 
tal j = state and 'ferromagnetic' interactions favoring fu- 
sion into a total j = I state. 

Unitary models- For the unitary case, such su(2)i; gener- 
alizations of the golden chain model were first discussed in 
Ref. |23] and further explored in Ref. [34] Like their su(2)3 
counterpart these generalized anyon chains are found to be 
exactly solvable and their gapless collective states can be de- 
scribed in terms of conformal field theories. To be specific, the 
critical theory for the generalized antiferromagnetic chains is 
known to be the minimal model A4{k + 1, k + 2) with central 
charge c = 1 — 6/((fc+ l){k + 2)), while for ferromagnetic 
chains it is known to be the Zk parafermion theory with cen- 
tral charge c = 2{k — l)/{k + 2). In both cases, the terms in 
the Hamiltonian are found to form a Temperley-Lieb algebra 
with isotopy parameter d — 2 cos(7r/ {k + 2)). 

Non-unitary models.- In turning to the Galois conjugated 
models of these unitary su(2)fc anyon chains, we first recall 
that for a given su(2)fc anyonic theory there are several allowed 
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TABLE ni. Overview of the gapless theories for the various chain models. 



q deformation parameters (see also Sec. II B i 



^ ^ gP-27rj/(fc+2) 



(24) 



as illustrated in Fig. |9] For each integer index 1 < p < 
{k + 2)/ 2 (and the further restriction that p and k + 2 are rela- 
tive prime) the corresponding q deformation parameter estab- 
lishes a solution of the pentagon equations of the underlying 
su(2)fc theory and defines the explicit form of the F-symbols. 
The unitary su(2)fe models discussed above all correspond to 
the first primitive root, i.e. p = 1, in this sequence. The 
other primitive roots (with p > 2) then give rise to the Ga- 
lois conjugated models with non-unitary i^-symbols and non- 
Hermitian microscopic Hamiltonians constructed using these 
i^-symbols. In particular, we note that for k ~ 2 there is 
no such Galois conjugated theory, for /c = 3 there is exactly 
one Galois conjugated theory (the Yang-Lee theory discussed 
above), and for fc > 4 there are in general multiple Galois 
conjugated theories. 

Like their unitary counterparts all Galois conjugated su(2)fe 
chain models permit an analytical solution. Since the line 
of argument closely follows the one given in Section [nC] for 
fc = 3, we will be brief here, and merely state the main argu- 



47ri/ll 



q = e 




27ri/ll 



► Re g 



ments and results. In order to simplify the presentation we will 
further focus our arguments to anyonic theories with odd fc 
which allows to restrict the discussion to those which possess 
only integer (not half-integer) values of "angular momenta" j. 
(Even values of fc work in an analgous manner.) The Galois 
conjugated su(2)fc chain models can then be directly mapped 
onto the height representations of certain restricted solid on 
solid models where the heights take the values 0, 1/2, . . . , fc/2 
in the Dynkin diagram Ak illustrated in Fig. [5] The Hamilto- 
nian terms acting in this height representation then again form 
a Temperley-Lieb algebra"*^ with an isotopy parameter which 
now turns out to depend on the sign of the interactions. For 
antiferromagnetic coupling this generalized isotopy parameter 
is found to be d = 2 cos(p7r/ (fc + 2)) and the conformal the- 
ory describing gapless interacting system is the non-unitary 
minimal model Ai{k + 2 — p, k + 2) with central charge^l 

For ferromagnetic couplings the iso- 
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(fc+2)(fe+2-p) 



topy parameter becomes d — 2 cos ^ ^*^"fc^_^2^''^ ) '■^^ 8^P" 
less theory is the non-unitary minimal model A4(p, fc+2) with 



central charge c — 1 



6{k+2~py 



In the case of the Galois- 



p(fe+2) 

conjugated models, the central charges are related to the quan- 
tum dimensions in the following way''^ ^* Writing the quan- 
tum dimension as d = 2cos(p7r/(fc + 2)) — 2cos{tt/5) 
(where (5 is a fraction for the Galois -conjugated models), one 
can obtain the central charge via c ~ 1 — s{S^-i) • Note that 
this is also the correct result for the anti-ferromagnetic unitary 
chain with p = 1, but the central charge for the ferromagnetic 
unitary chains is given by c = 2 — 6/(5. 



FIG. 9. The g-deformation parameters for a general su(2)fc theory 
correspond to different primitive roots of unity. Only the first primi- 
tive root in this sequence, e.g. q = e^'^'''^^ in the illustrated example 
of fc = 9, gives rise to a unitary model. 



Our results for unitary and non-unitary models are summa- 
rized in Table III We note that for p = 2 and fc = 3 we exactly 



reproduce the results of the Yang-Lee chains discussed above. 
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m. DOUBLED YANG-LEE MODELS 

A. The ladder Hamiltonian 

In this second part of the manuscript, we turn to 'quan- 
tum double' variants of the anyonic chains discussed in the 
first part. The unitary incarnations of these quantum double 
models have been introduced in the context of exotic quan- 
tum phase transitions in time-reversal invariant systems that 
are driven by topology fluctuations^^. The specific model is 



defined on a ladder geometry, shown in figure 10 
The Hamiltonian 



-ff ladder 



l{r),l 



rungs r 



plaquettes p 



(25) 



consists of two competing terms. 

The first term favors the trivial label 1 on each rung of the 
ladder, while the second term favors the no-flux state for all 
plaquettes. As shown in Ref. [28] the projector onto the flux 
through a square plaquette can be expressed in terms of the 
unitary/non-unitary F-matrices (|6|/(|7]i. This term is equivalent 
to the plaquette term in the Levin-Wen modelJ^, which are 
defined on a different lattice, the honeycomb lattice. 

Explicitly, the plaquette term reads 



30(p),i 
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(26) 



where dg denotes the quantum dimension of particle type s, 
i.e. di = 1 and (1^- ~ cf). D denotes the total quantum di- 
mension, D — d\ + d'^ = \J2 + (/) for Fibonacci (as well 
as Yang-Lee) anyons. The latin and greek labels denote the 




FIG. 10. The geometry of the Fibonacci and Yang-Lee ladders. To 
discuss the physics, it is enlightening to 'thicken' the ladder, and 
consider the two dimensional surface thus obtained. 



solvable point 
(c = 8/35) 




solvable point 
(c = -44/5) 



FIG. 11. Phase diagram of the doubled Yang-Lee chain. 



degrees of freedom, and any of these takes one of the values 
{1,t}. We note that the Hilbert space of the ladder mod- 
els consists of all possible labelings of the rungs and the legs, 
such that at each vertex, the Fibonacci fusion rules are obeyed. 

With this description of the ladder models, we can easily go 
back and forth between the Fibonacci anyon ladder, and the 
Yang-Lee anyon ladder, simply by choosing the corresonding 
set of i^-symbols, namely equations (|6]l or (|7]) respectively. 



B. The phase diagram 



To discuss the phase diagram of both the original and Ga- 



lois conjugated model, shown in figure 1 1 we parametrize the 
couplings Jr — sin d and Jp — cos B in terms of an angle Q. 

We start with the first gapped phase, for 7r/4 < < tt, most 
easily discussed at the special point 6* = 7r/2, or J,. = 1 and 
Jp = 0. In the ground state at this point, all rungs have the 
trivial label. In the two-dimensional surface geometry, this 
means that no r-fluxes go through the rungs, implying that 
the rungs can be completely pinched off, changing the geom- 
etry to that of two, disconnected cylinders, one for each leg of 
the ladder. Each of the cylinders accommodates two ground 
states, either with or without r-flux, leading to a four-fold 
ground state degeneracy. The lowest (gapped) excited state 
consists of configurations where one rung of the ladder 'con- 
tains' a T flux. 

In the other extended gapped phase, for — 7r/2 < ^^ < 7r/4, 
let us discuss the special point ^? = 0, or Jp = 1 and J,. = 0. 
Here the situation is reversed, and no x-fluxes go through the 
plaquettes. Thus they can be closed of, giving rise to a geome- 
try consisting of a single cylinder Again, this cylinder can ac- 
commodate two ground states, with or without flux, resulting 
in a two-fold degenerate ground state. The lowest (gapped) 



11 



excited state consists of configurations with a single T-flux 
going through a plaquette, effectively piercing a hole through 
the cylinder. 



Precisely at the point where both couplings are equal in 
strength, the gap closes, and the system is critical. At this 
point, the geometry is fluctuating at all length scales, interpo- 
lating between the two extremes of having one or two cylin- 
ders, respectively. In addition, also precisely at this point, the 
(critical) model is exactly solvable, as explained in the next 
section. 



Finally, for tt < < 37r/2, there is an extended critical re- 
gion, which is characterized by another exactly solvable point, 
atO = 57r/4, where both couplings are again of equal strength, 
but negative. 



C. Analytical solution 

The analytic solution of the Fibonacci ladder model at the 
two exactly solvable points has been described in detail in (the 
supplementary material of) reference |28] Thus, we will be 
rather brief here, and point out some of the crucial steps of the 
mapping of the ladder model onto an exactly solved (restricted 
solid-on-solid) model. We will then quickly mention in which 
way this solution of the unitary Fibonacci case is changed, if 
one takes the Galois conjugate, to go to the case of the Yang- 
Lee ladder model. The most important step in the analytic 
solution of the model, is to perform a basis transformation on 
the Hilbert space. After this transformation, the model can 
be mapped onto a restricted solid-on-solid model, where in 
this case, the allowed height variables lie on the vertices of 
the Dynkin diagram of the Lie algebra Dq, as it will turn out. 
This model is exactly solvable, both in the original, as well as 
the Galois conjugated case. 



The basis transformation: At each rung, one applies an F-transformation, changing the basis building blocks of the ladder: 

fli 0,3 ai 03 05 



C4 



hi 



bi 



FIG. 12. Basis transformation of the ladder model. 



Thus, we can express the Hilbert space in terms of the vari- 
ables {a2i-i, 621-1, C2i}, or, equivalently, in terms of the vari- 
ables {a2i-i, b2i-i, d2i}, as displayed in figure 12 The map- 



ping to the restricted solid-on-solid model is performed in 
terms of the latter variables. In particular, one can determine 
which 'values' the links can take on the even and odd sublat- 
tices. The variables on the even sublattice d2i can take the 
values while the combined variables (a2i-i, 622-1) 

can take the values {(1, 1), (1, r), (t, 1), (r, x)}, but there 
are constraints on which values of (a2i-i, 52i_i) and d2i can 
occur next to each other. From the fusion rules, it easily fol- 
lows that the combination of variables which are consistent 
with one another, have to be adjacent on the 'height graph', as 



shown in Fig. 13 




(1,1) 



(l,r) 



(t,1) 



FIG. 13. The allowed 'height' configurations of the ladder models. 
This diagram is the Dynkin diagram of the Lie algebra — so( 12). 



From now on, we will call the allowed variables (or pairs 
of variables on the odd sublattice), heights. The Hilbert space 
of the ladder model consists of all height configurations on a 
chain of length L, in such a way that two heights on neigh- 
boring sites also have to be neighbors (i.e., be connected by 
a solid line) in the height diagram of Fig. [T3] In this way, 
the model can be mapped on the Dg version of the restricted 
solid on solid model, and precisely when Jr — Jp ~ ±1, 
the model is critical, and can be solved exactly. We will not 
go into the details here (which can be found in Ref. 28 ), but 
we will quickly state the result for this critical behavior in 
the original, Fibonacci anyon setting, followed by the results 
on the Galois conjugated Yang-Lee ladder The details of the 
phase diagram, which closely mimics the phase diagram of 
the unitary ('not Galois-conjugated') ladder, and the numeri- 
cal results will be presented in the next subsection. 

As was the case for the chain models, also in the ladder case 
do the parts of the Hamiltonian form a representation of the 
Temperley-Lieb algebra, but in this case, the algebra is tha^ 
based on the Dynkin diagram D^. The possible associated val- 
ues of the rf-isotopy parameter are given by d = 2 cos(^^^), 
where ft, = 10, and the nij are the exponents of the Lie algebra 
Dq, namely rrij = 1, 3, 5, 7, 9, 5. In the original ladder model 
nij = 1, for both cases Jr — Jp — ±1- The associated value 
of the d-isotopy reads d = 2 cos(7r/10) = ^/2 + cj). In the an- 
tiferromagnetic case, i.e. or in the parametriza- 
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1/40 247/280 1287/280 


2 


-11/40 




Yang-Lee ladder 




3 


2/5 9/35 104/35 


3 


-2/5 


coupling 


CFT 


central charge 


d-isotopy 


4 


9/8 -1/56 95/56 


4 


-3/8 


AFM (e = 7r/4) 


M{7,IQ) 


c = 8/35 


2 cos (If) 


5 


11/5 2/35 27/35 


5 


-1/5 


FM {e = 57r/4) 


M(3,10) 


c = -44/5 


2cos(f^) 


6 


29/8 27/56 11/56 


6 


1/8 
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27/5 44/35 -1/35 


7 


3/5 


TABLE IV. Gapless theories for the ladder models. 
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tion Jp = cos 



Jr = sin 6, the point 9 = 7r/4, the critical 
behavior is described by the minimal model Ai{9, 10), which 
has central charge c — 14/15. With the tilde, we d enote 
that we are considering the (AgjDg) modular invarianiSEZl 
instead of the 'usual' (diagonal) {As,Aq) invariant (we note 
that the (Ag, Dq) invariant minimal model (9, 10) is equiv- 
alent to the coset model (62)1 x (G'2)i/(G'2)2), where {G2)k 
denotes the Wess-Zumino-Witten CFT at level k based on the 
Lie-algebra G2. 

We will now switch our attention to the Galois conjugated 
version of the ladder model. The phase diagram of the Yang- 
Lee ladder is depicted in figure 1 1 and the location of the 
gapped phases, as well as the critical behavior is exactly the 
same as in the Fibonacci ladder model. Therefore, we will 
be somewhat brief in the description of these phases (see 
the introduction of this section for a brief description of the 
gapped phases), and focus on the differences, which occur at 
the exactly solvable points, namely 9 = 7r/4, 57r/4, where 
the model is critical, but with a different critical behavior in 
comparison to the Fibonacci ladder. 

The anti ferromagnetic {9 — 7r/4) critical point.- As was 
the case in the Yang-Lee chains, the critical behavior of the 
Galois conjugated model is obtained by picking a different 
d-isotopy parameter, corresponding to the Galois conjugate 
under consideration. For Fibonacci anyons, there is only 
one Galois conjugate, but nevertheless, there seem to be dif- 
ferent possible values for the d-isotopy parameter, namely 
d = 2 cos(^^^), with h ~ 10, and raj — 1, 3, 5, 7, 9, 5. In 
fact, the value we need is nij — 3, which corresponds to the 
value obtained from the Fibonacci case, under Galois conju- 
gation, namely d — ■\/2 — 1/0 = 2 cos(37r/10). 

It is a simple matter to check that this value corresponds to 
a model with central charge c — 8/35, by using the same ar- 
guments as we gave in section [TTfI This central charge corre- 
sponds to the central charge of the non-unitary minimal model 
A^(7, 10). Just as in the unitary casJ^ll the critical CFT de- 
scribing the Yang-Lee model at 6* = 7r/4, is not the usual ('di- 
agonal') modular invariant, but rather the (A, D) invariant, or, 
to be more precise, {Aq,Dq). Sed^for information on the 
different modular invariants of the minimal models. We will 
denote this model by A^(7, 10). The field content of these 
models is given in tables [V] and | VI| below. Details of the spec- 
trum of this model will be discussed in the section llll DI below. 



TABLE V. Kac table for the non-unitary minimal models Al(7, 10) 
and A1(3, 10). Only odd values of r are displayed, so that each field 
appears only once. 



The ferromagnetic [9 — b-K /A^) critical point.- In the non- 
unitary case, we can obtain the critical theory for the opposite 
sign of the interaction, by simply swapping the sign of the 
d-isotopy parameter, which hence corresponds to the value 
d — 2cos(77r/10) = —\J2— l/4>, from which it is a sim- 
ple matter to extract the central charge, which is given by 
c ~ —44/5, corresponding to the minimal model A^(3, 10). 
This minimal model also comes in two incarnations, one cor- 
responding to the modular invariant (742,^9), the other to 
{A2tDq). (The field content of these models is also given 
in tables [V] and VI below.) The latter is realized in the Yang- 
Lee ladder at 6* = 57r/4, and we will describe the spectrum 
in more detail below. Before we do that, however, we will 
first summarize the critical theories of both the Fibonacci and 



Yang-Lee anyon ladders in table IV 



M(7, 10) 



>1(3, 10) 



h{r, s) 


s = 1 3 5 


/i(r, s) 


s = 1 


r- = 1 


13/7 4/7 


r = 1 





3 


2/5 9/35 -1/35 


3 


-2/5 


5 


11/5 2/35 27/35 


5 


-1/5 



TABLE VI. Kac table for the '{A, Z))-modular invariant' non-unitary 
minimal models M{7, 10) and A1(3, 10). The fields with the con- 
formal dimensions in bold (i.e. those with label s = 5) appear twice. 



D. Numerical results 

We finally present the numerical spectra of the conjugated 
ladder model at the two critical points discussed in the previ- 
ous section. The spectrum for the critical point ?a9 — 7r/4, is 
given in Fig. [14} where we indicated the locations of the pri- 
mary fields of A^(7, 10) (given in table VI 1 by green squares, 
as well as some low-lying descendants with red circles. As 
usual, there are only two free parameters to match the numer- 
ically obtained spectrum with the result obtained from con- 
formal field theory, so the fact that the six lowest primaries, as 
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FIG. 14. Conformal energy spectrum of the critical points in the 
doubled 'antiferromagnetic' Yang-Lee chain. The spectrum matches 
the non-unitary minimal model M (7, 10) with central charge c = 
8/35. Primary fields of the conformal field theory are indicated by 
green squares, descendant fields by red circles. 
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FIG. 15. Conformal energy spectrum of the critical points in the 
doubled 'ferromagnetic' Yang-Lee chain. The spectrum matches 
the non-unitary minimal model M(3, 10) with central charge c = 
—44/5. Primary fields of the conformal field theory are indicated by 
green squares, descendant fields by red circles. 



well as several descendants match to high precision (limited 
by finite size effects) is a very non-trivial check on our results. 



In Figure 15 we give the numerical spectrum of the ferro- 



magnetic Yang-Lee ladder, at 6* = 57r/4, which is character- 
istic of the critical phase extending over 9 G {tt,3tt/2). In 
this case the critical behavior is described by the A4{3, 10) 
non-unitary conformal field theory, and, as for the antiferro- 
magnetic case, we were able to identify the primary fields, as 
well as several low-lying descendant fields, as indicated in the 



IV. DISCUSSION & SUMMARY 

In this paper, we studied the collective states of Yang-Lee 
anyons, a family of non-unitary, non-Abelian anyons which 
are close cousins of the unitary Fibonacci anyons. Non- 
unitary anyons of this form have attracted interest in the con- 
text of studies of certain quantum Hall wavefunctions, includ- 
ing the Gaffnian state Both Yang-Lee and Fibonacci anyons 
arise from the same anyonic theory, su(2)3, and in particular 
they share the same fusion rules. The key distinction between 
the two anyon types is that Yang-Lee anyons are non-unitary 
and relate to their unitary counterparts, the Fibonacci anyons, 
via 'Galois conjugation'. We have generalized this concept to 
arbitrary su(2)fc anyonic theories. 

To characterize the collective states formed by a set of 
anyons in the presence of pairwise interactions, we have 
considered one-dimensional models of interacting Yang-Lee 
anyons similar to the golden chain model of the unitary case^^] 
Analogous to the case of interacting Fibonacci anyons, the 
collective states of such chains of Yang-Lee anyons are found 
to be critical and the gapless theories are described by cer- 
tain non-unitary conformal field theories (which depend on 



figure. The fields of the CFT are given in table VI 



the sign of the coupUng, see Table III 1. 

The non-unitary chain models are found to exhibit rather 
peculiar features related to the presence of a non-local, topo- 
logical symmetry first observed in the unitary models. For 
the multicritical gapless theories of these chain models this 
topological symmetry, whose related symmetry operator com- 
mutes with the Hamiltonian, allows to classify all operators by 
a topology symmetry sector. It turns out that all relevant oper- 
ators corresponding to (uniform) perturbations of the gapless 
system are in a different symmetry sector from the one of the 
ground state. In the unitary case, this mechanism which ef- 
fectively protects the gapless ground state from local pertur- 
bations has led to an interpretation^'* of these gapless modes 
of a chain of interacting anyons as edge states of the par- 
ent topological liquid of which the anyons are excitations of. 
This has further led to the conclusion that interactions between 
anyons (in two-dimensional arrangements) result in a splitting 
of the macroscopic degeneracy of a set of non-Abelian anyons 
and the nucleation of a new topological liquid inside and spa- 
tially separated from the parent liquid of which the anyons are 
excitation^^^l^lJ 

For the non-unitary models, studied in this manuscript, a 
crucial distinction comes in the symmetry sector assigned to 
the ground state. While in the unitary case, the ground state 
was found to be in the flux-free (topological trivial) sector, we 
find that the ground state in the non-unitary case exhibits a 
non-trivial topological flux (i.e. is in the non-trivial topolog- 
ical symmetry sector). It remains an open question whether 
such a spontaneous creation of topological flux can occur in 
the ground state of a system of interacting anyonic quasipar- 
ticles residing in a bulk-gapped topological liquid, which has 
no topological flux associated with it. This may be an indica- 
tion that the non-unitary anyons are not massive quasiparticles 
of a gapped quantum liquid, but in fact excitations of a gapless 
quantum liquid. 

In summary, we have investigated one-dimensional models 
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of non-unitary anyons based on su{2)k anyonic theories. The 
collective ground states are found to display critical behavior, 
similar to their unitary counterparts to which they are related 
via a 'Galois conjugation'. We described in detail the non- 
unitary conformal field theories capturing their critical behav- 
ior and commented on possible physical implications of our 
results in connection with proposed quantum Hall states re- 
lated to such non-unitary CFTs. In addition, we also studied 
the phase diagram of an interacting ladder model, which is 
a quasi one-dimensional version of the Levin- Wen model for 
non-unitary Yang-Lee anyons. 
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Appendix A: Detailed description of the conformal energy 
spectra. 

In this section, we will describe in some detail the structure 
of the descendant fields present in the numerically obtained 
spectra. In the case of the Yang-Lee ladders, this structure 
is somewhat different from the 'usual' structure. The object 
containing the information about the number of states at a spe- 
cific energy and momentum (in the thermodynamic limit) is 
the partition function, which we will denote by Ztot ■ 

Modular invariance of the partition function on the torus 
constraints the possible partition functions of rational confor- 
mal field theories at a particular central chargd^. In this pa- 
per, we will only be concerned with the partition functions of 
minimal models, which are described in terms of the chiral 
characters associated to the primary fields. An explicit ex- 
pression for these chiral characters associated to the primary 
fields of the minimal models M{p,p') is given in, for 

instance, chapter 8 of reference l60l 



„^,(P,P')^„^ _ g \^ ( Hkpp'+rp-sp') _ (kp+s)(kp'+r) 



where (g)oo = nfe>o(l ^ ^"^^ ^(r.s) 



the conformal dimension of the field 



(Al) 

_ jpr-p' sf-{p-p'y 
App' 



Dropping the labels {p,p'), all the possible partition func- 
tions can be written as 

Ztot = ^ ^'I(r,s);(r' ,s')C\r,s)C^lr' ,s') ■ (A2) 
(r,s);(r',s') 

For the minimal models, there is always the so-called 
'diagonal modular invariant', for which M is diagonal. 



Af, 



(r,s);(r',s') 



5r.r'5s,s'^ and the sum in the partition func- 



exist other modular invariant partition functions. The criti- 
cal points of the ladder models realize some of these, as we 
pointed out in the main text. 

In terms of the Virasoro generators Lq and Zo, the Hamilto- 
nian and momentum operators can be written (after an appro- 
priate rescaling and shift) as H = Lq+Lq and P = (Lq—Lq). 
The partition function can be written as a trace over the 
Hilbertspace, Ztot — tr q^^q^", which allows us to extract 
both the energies and momenta for each state. 

In particular, the energy of a state is given by the sum of 
the powers of q and q. We note that this energy is measured 
with respect to the energy corresponding to the energy of the 
identity operator 1, which we have set to zero. In addition, 
the momenta of a state is given by Lq — Lq, in units of 2ti/L, 
and measured with respect to the momenta of the correspond- 
ing primary field, which is not determined by conformal field 
theory, but rather by the specific Hamiltonian realization. 



1. The ferromagnetic Yang-Lee ladder 

The critical theory describing the ferromagnetic Yang-Lee 
ladder is the (Dg: ^2) modular invariant A1(3, 10). The total 
partition function of this theory can be written as followJ^ 
(we will drop the labels {p,p') — (3, 10) from the characters 
for convenience) 

Z = |ch(i,i) |2 + |ch(3,i) |2 + 2|ch(5,i) P + |ch(7,i) P (A3) 

+ |ch(g_i)p + (ch(i^i)ch(i 2) + ch(3,i)ch(3 2) + c.c.) . 

By expanding the partition function, in terms of q and q, 
we can completely explain the low-lying part of the spectrum 



displayed in figure 15 Let us start with the vacuum sector 
of the theory, which corresponds to the part of the partition 
function with integer powers of q and q: 

Zo = |ch(i,i) P + |ch(9_i) p + (ch(i_i)ch*i 2) + c-C-) = 
1 + 2g2 + 2g2 + . . . . (A4) 

The total power of q and q gives the energy of the states, 
while the difference between the powers gives the momentum. 
Thus, the character of the vacuum sector implies the presence 
of two states at energy E — 2 and = 2 as well as two 
states at energy E ~ 2 and momentum = —2. Note that 
we count the states irrespective of their Ky sector. 

We will continue with the descendants corresponding to the 
primary field with energy E — —2/5 and momentum ~ 0. 
At this energy and momentum, there are in fact two states. 
This sector of the partition function reads 



^-2/5 



2|ch(5,i)P 



2q-^/^q-^/^ (l + q + q + 2q^+ 2f 



■) , (A5) 



which implies two states at {E, K^) = (1 — 2/5, 1) and two at 
{E,K^) = (1-2/5,-1), as well as two states at (2 -2/5,0) 
and four at both (2 - 2/5, 2) and (2 - 2/5, -2), all of which 
is reproduced (modulo finite size effects) in the spectrum in 



tion runs over all primary fields. However, in general, there figure 15 
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Finally, we consider the descendants corresponding to the 
primary field at E = —4/5, which also has momentum = 
0. This sector of the partition function reads 



Z_4/5 = |ch(3^i)p + |ch(7_i)p + (ch(3^i)Ch(3 2) + C.C.) 



^-2/5 --2/5 fi_^2q + 2q + 3q^ + Sq^ + Aqq 



(A6) 



This implies the presence of two states at both {E,Kx) = 
(1 - 4/5, 1) and (1 - 4/5, 1), as well as three states at both 
(2 - 4/5, 2) and (2 - 4/5, -2). hi addition, there should be 
four states at (2 — 4/5, 0). This seems to be at odds with the 



figure 15 but it turns out that the state at = and Ky ~ tt, 
and energy of approximately £' w 2 — 4/5 (or, more precisely, 
E w 1.2144495) is in fact doubly degenerate, so indeed, there 
are four states, as expected from the partition function. 



2. The antiferromagnetic Yang-Lee ladder 

A similar analysis can be performed for the antiferromag- 
netic critical point of the Yang-Lee ladder. We will be brief 
here, because the analysis is identical to the one for the ferro- 
magnetic critical point. 

There are nine sectors in this particular case, corresponding 



to the fields displayed in table VI so the partition function can 
be written as 



Zs/7 



^4/5 ■ 



z. 
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8/35 
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2/35 



^54/35 



(A7) 



Although we will not consider all of these sectors (some give 
rise to states which are too high in energy to be unambigu- 
ously identified, due to finite size effects), we will neverthe- 
less give the form of the partition function restricted to each of 
these sectors. We will drop the labels {p,p') — (7, 10) from 



the characters 
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The parts of the partition function read explicitly 



Za = l + q^ + q^ + 
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^54/35 = 
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q + iq^ + qq + W + . . 

2 I /I„;t I 1^2 



^9/35-9/35 (1 + 2q + 2g + 3(7^ + Aqq + 3^ 

2^2/35-2/35 A 



q + q + 2q^ + qq + 2f + ..) 
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2^27/35-27/35 (^i + q + q + 2q^ + qq + 2q^ + . . .) 



With this information, it is rather straight forward to check 
that the primaries and descendant fields we indicated in figure 
14 indeed come with the right multiplicities. 



3. The ferromagnetic Yang-Lee chain 

After having dealt with the (anti)ferromagnetic Yang-Lee 
ladders in quite some detail, we will content ourselves here by 
giving the (diagonal) partition functions describing the spec- 
tra, and note that for all states where finite size effects allow 
us to make an identification, we obtain full agreement with the 
CFT prediction. 

The partition function in case of the ferromagnetic Yang- 
Lee chain, the critical model is given by the diagonal invariant 
of the model Ai{2, 5), whose partition function reads (drop- 
ping the label (2, 5) on the chiral character) 



Z = |ch(i,i)|2 + |ch(3,i)|2 = Zo + Z_2/5 



(A8) 



with 



Zo = l + q"^ + q^ + q^ + f + q"^ + q^f + q^+ 
q^ + q^f + q^f + f+ 
2q^ + q^q^ + q^q^ + q^q"" + 2f + ■■■ 



(A9) 



^-2/5 - q'^'^^r^'^ {l + q + q + q^ +qq + q^+ (AlO) 



q"^ + q^q + qf +(f~ 



2q* + q^q + q^f + qf + 2^+ 

2q^ + 2q'^q + q^f + q^q^ + 2qt + 2^ + 

3q^ + 2q^q + 2q^q^ + q^q^ + 2gV + '2qf + 3q^ 



All these states, and in fact quite a few more which we did not 
give here, are reproduced in the spectrum of the ferromagnetic 
Yang-Lee chain, as shown in figure [8] 



4. The antiferromagnetic Yang-Lee chain 

Finally, we deal with the antiferromagnetic Yang-Lee 
chain, the critical behavior of which is described by the (diag- 
onal invariant of the) minimal model A4{3, 5), which has four 
primary fields, as described in the main text. The relevant par- 
tition function is given by (again, dropping the labels (3, 5) on 
the chiral characters) 



|ch(i 



|c]i(. 



3,1)1 



|ch( 



1,2)1 



|cll(3,2)P 
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= Zq + Z2/5 + ^3/2 + ^■^i/io , (Al 1) + 2(z3 + g2 - + ^ . . . ) 

with 



""^'^^''^ ^ ^ We should note that, as can be seen in figure]?] these fields are 

Z2/5 = q^/^g^/^ (1 + q + ^ + 2q^ + qq + 2q^ + • • • ) reproduced, and that (for L — 2,2 sites), the fields with confor- 

_3/4 3/4/, - \ mal dimension ft, = —1/20 and /i = 3/4 occur at momentum 

3/2 ==9 9 (1 + 9 + g + • • • j /^^ = 0, while the identity fields and the field with conformal 

^-1/10 = + q + q + q^ + qq + q'^+ dimension h = 1/5 occur at momentum = tt. 
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